Molecular Spiders in One Dimension 
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Molecular spiders are synthetic bio-molecular systems which have "legs" made of short single- 
stranded segments of DNA. Spiders move on a surface covered with single-stranded DNA segments 
complementary to legs. Different mappings are established between various models of spiders and 
simple exclusion processes. For spiders with simple gait and varying number of legs we compute the 
diffusion coefficient; when the hopping is biased we also compute their velocity. 



I. INTRODUCTION 

In recent years, chemists have constructed a number of 
synthetic molecular systems which can move on surfaces 
and tracks (see e.g. [H, [3, [1] and a comprehensive review 
One such object is a multi-pedal molecular spider 
whose legs are short single-stranded segments of DNA 
These spiders can move on a surface covered with 
single-stranded DNA segments, called substrates. The 
substrate DNA is complementary to the leg DNA. The 
motion proceeds as legs bind to the surface DNA through 
the Watson-Crick mechanism, then dissociate, then re- 
bind again, etc. More precisely, a bond on the substrate 
with an attached leg is first cleaved [3], and the leg then 
dissociates from the affected substrate (which we shall 
call product). The leg then rebinds again to the new 
substrate or to the product leading to the motion of the 
spider. 

The rate of attachment of a leg of the spider to the 
substrate and the rate of detachment from the substrate 
are different from the corresponding rates involving the 
product instead of the substrate. Hence for the proper 
description of the motion of a single spider one must keep 
track of its entire trajectory. This memory requirement 
makes the problem non-Markovian and generally 

intractable analytically even in the case of a single spider. 
(We shall address this problem in a separate paper [1].) 
Many interacting spiders add another level of complexity. 
Even if the rates were the same for the substrate and the 
product, the properties of the spider (e.g., the nature 
of its gait) represent another major challenge. Here we 
separate this latter issue from the rest: we investigate 
how the structure of the spider affects its characteristics 
(velocity and diffusion coefficient). Further, we consider 
spiders with idealized gait — the goal is not to mimic 
complicated (and poorly known) gait of molecular spiders 
but to show that spiders' macroscopic characteristics are 
very sensitive to their structure and gait. 

We shall mostly focus on a single spider moving on a 
lattice. We shall also assume that the rate of attachment 
greatly exceeds the rate of detachment. In this situation 
the relative time when one leg is detached is negligible 
and hence the possibility that two or more legs are de- 
tached simultaneously can be disregarded. Therefore the 
present model posits that spiders remain fully attached 



and never leave the surface % (This differs fr om the 
actual situation when a few legs may be simultaneously 
detached.) 

The spiders are defined as follows. Legs can jump in- 
dependently at constant rates if they do not violate the 
restrictions below. We mainly consider symmetric spi- 
ders where we set all these rates to one, or biased spiders 
whose legs can only move to the right at rate one, but 
some special gaits are also investigated. The fundamental 
restriction on the spider's motion is the exclusion princi- 
ple: Two legs cannot bind to the same site. Additional 
constraints keep the legs close to each other. We mainly 
consider two types of spiders with the simplest feasible 
constraints: 

Centipedes (or local spiders). A leg of a centipede 
can step to nearest neighbor sites provided that it re- 
mains within distance s from the adjacent legs. (This 
threshold is assumed to be the same for each pair of ad- 
jacent legs.) 

Spiders (or global spiders). Legs of these spiders can 
step to nearest neighbor sites as long as all legs remain 
within distance S. 

The above properties of the gait guarantee that in one 
dimension the order of the legs never changes. The above 
constraints seem equally natural in one dimension, while 
in two dimensions the global constraint appears more rea- 
sonable. 

We shall also briefly discuss a third type of spider 
where the nearest neighbor restriction on the hopping is 
relaxed. For these quick spiders, legs can step anywhere 
within distance S from all legs. Quick spiders have been 
proposed and studied numerically in Ref. p^ . 

The above assumptions about the gait and the dis- 
regard of memory effects leave little hope for quantita- 
tive modeling, but simplicity can help to shed light on 
qualitative behaviors. Therefore we study in depth a 
single spider with aforementioned gait moving on a one- 
dimensional lattice, and more briefly probe the influence 
of the gait and many-spider effects. 

The rest of this paper is organized as follows. In Sec. [Ill 
we analyze bipedal spiders (i.e., spiders with two legs). 
This framework provides a useful laboratory to probe 
various techniques. Bipedal spiders also closely resem- 
ble molecular motors and the methods developed for 
studying molecular motors are fruitful for studying indi- 
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vidual spiders Sec. lIIII we examine multi-pedal 

spiders (i.e., spiders with L > 3 legs). We show that the 
spider with local constraint and s = 2 is isomorphic to a 
simple exclusion process (SEP) on a line with L — 1 sites 
and open boundary conditions; an even simpler isomor- 
phism exists between spiders with global constraint and 
the SEP on the ring. These connections allow us to ex- 
tract some spider characteristics from results about the 
SEP. Quick spiders are briefly investigated in Sec. [iVl In 
Sec. |V] we show that the behavior of many interacting 
spiders can also be understood, at least in the practi- 
cally important situation of low spider density, via the 
connection with the SEP. Finally we stress limitations of 
our analysis and discuss possible extensions in Sec. IVII 



II. BIPEDAL SPIDER 

For bipedal spiders, the local and global constraints 
are equivalent, s = S. For the simplest mobile bipedal 
spider, the allowed distance between the legs is one or two 
lattice spacing, i.e., s = 2. Two possible configurations 
are (up to translation) 



and 



o • o • o 



(1) 



where we denote empty sites by 'o' and filled sites (to 
which the legs are attached) by There are obvious 
back and forth transitions between these configurations: 



and 



For symmetric spiders each leg jumps at rate one when 
possible, hence all the above four elementary moves hap- 
pen at rate unity. The diffusion coefficient of this bipedal 
spider is 



(2) 



To put this in perspective, we note that the diffusion 
coefficient of a random walker which hops to the right 
and left with unit rates is D = 1. Thus adding a leg and 
requiring the legs to stay within distance two to each 
other reduces the diffusion coefhcient by a factor 4. 

Generally for symmetric bipedal spiders with arbitrary 
s, there are s possible configurations labeled by the 
inter- leg distance, £ — 1, . . . , s. The transitions are 

Ci ^ 62 ^ . . . ^ 

The diffusion coefficient Dg of this bipedal spider is 



1 / 1 

Ds = - [1- - 

2 V s 



(3) 



The above results apply to symmetric bipedal spiders 
which hop to the left and right with equal rates. Molec- 
ular motors usually undergo directed motion pT| . and 
one of the goals of future research is to control spiders to 



move preferentially in a certain direction. Here we ana- 
lyze such directed motion theoretically. For concreteness, 
we focus on the extreme bias when each leg can only hop 
to right at rate one. For instance, for the bipedal spider 
the most compact configuration evolves via ••o •o» ; 
the process o • • • o • involves hopping to the left 
and therefore it is forbidden in the biased case. 

For biased bipedal spiders the velocity and the diffu- 
sion coefficient are given by 



(4) 



In this section we give a pedestrian derivation of ([3]). 
The expressions Q for velocity and diffusion coefficient 
can be derived by utilizing the same technique; instead, 
we shall extract them from more general results for lame 
spiders fSec. lIICp . 

To set the notation and to explain how we compute the 
diffusion coefficient we begin with a random walk (which 
is a one-leg spider). Let P„(i) be the probability that 
the random walker is at site n at time t. This quantity 
evolves according to 



dPn 
dt 



Pn-l + Pn 



2Pr, 



(5) 



One can solve this equation and then use that solution to 
extract the diffusion coefficient. In the case of the spiders, 
however, master equations generalizing ([5]) are much less 
tractable, and therefore a more direct way of computing 
the diffusion coefficient is preferable. Here we describe 
one such approach 7]. It involves two steps. First, one 
ought to determine the mean-square displacement 



E 



n'Pr,. 



n— — QQ 



Then the basic formula 7] 



D = hm 

t—>oo 



2t 



(6) 



(7) 



allows to extract the diffusion coefficient. 

For the random walk, the mean-square displacement 
evolves according to 

, oo 

n^Pn-i + Pn+i - 2Pn) (8) 

n— — OO 

Transforming the first two sums we obtain 

OQ OO 

Ti— — OO n— — OO 

These identities ahow us to recast ^ into 

, OO 

lix') = [{n+ir + in-ir-2n']P,, 

n— — oo 

OO 

= 2 ^ P„ = 2 
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where the last equahty follows from normalization. Thus 
{x^) — 2t. Plugging this into ([7]) we recover the diffusion 
coefficient of the random walker D ^ 1. 

We now turn to the bipedal spider. We shall examine 
in detail only symmetric hopping. 



A. Bipedal Spider with s — 2 

For the bipedal spider with s = 2 there are two possible 
spider configurations. Denote by Pn(t) and Qn{t) the 
probabilities that at time t the spider is in respective 
configurations ([T]), namely 



Pn = Prob[. .] , Q„ = Prob[. o .] , 



(10) 



with the left leg being at site n. The governing equations 
for these probabilities are 



dPn 
dt 

dQn 
dt 



(lla) 
(lib) 



The mean position of the legs or the 'center of mass' of 
the spider in a configuration corresponding to P„ (resp. 
Qn) is located at n + ^ (resp. n + 1). Thus the mean- 
square displacement is 



E 



n— — oo 



" + ^ ) Pn + in + if' 



(12) 



and it evolves according to 



dt 



J2 ("+^) (Qn + Qn-l - 'iPn 
oc 

f J2 {n+lf{Pn+l+Pn-2Qn) 



Utilizing the same tricks as in ^ we recast the above 
equation into 



dt 



OO 



(13) 



The last identity is implied by normalization and its va- 
lidity also follows from Eqs. ()llap - (|llbp . Integrating (flS]) 
yields (a;^) = ■^t which in conjunction with ^ leads to 
the previously announced result, Eq. 



General Case 



In the general case (s > 2) we denote 
Pn — Prob[»^-^»] 



(14) 



the probability to occupy sites n and n + £. These prob- 
abilities obey 



_ p2 , p2 _ 9 pi 
dt -P-^+P- 

dP' 



dt 
dP, 



(15a) 



ptt\ + Pt^\ + pt' + Pt'^^Pi (15b) 



n _ ps— 1 I ps-1 _ 9 ps 

dt -^n+l+^n ^^n 



(15c) 



where equations (|15b|l apply for 2 < ^ < s — 1. The 
mean-square displacement is given by 



E E 

n— — DO £—1 



n + - I PI 



(16) 



Using Eqs. p5ap - p5cp and applying the same tricks as 
above to simplify the sums, we obtain 



e=i 



(17) 



where we — J2n weight of configurations of the 

type . The sum on the right-hand side of Eq. (fT7|) is 
equal to one due to normalization. To determine wi and 
Ws one does not need to solve an infinite set of the master 
equations l|15ap - (|15cp . Instead, we take Eqs. (|15aP " (|15cp 
and sum each of them over all n to yield a closed system 
of equations for the weights 



dwi 

~dr 

dwe 
IT 
dws 
IT 



: 2(W2 - Wl) 

2{we-i + we+i — 2we) 

2{Ws-l - Ws) 



(18a) 
(18b) 
(18c) 



If initially wi = . . . = Ws = 1/s, then Eqs. (fT8a|l - (fT8cll 
show that this remains valid forever. Even if we start 
with an arbitrary initial condition, all the weights Ws 
relax exponentially fast toward the 'equilibrium' value 
1/s. Thus the right-hand side of ([T7]) becomes 1 — 1/s 
yielding (x^) = (1 — 1/s) t which in conjimction with ^ 
leads to Eq. ©. 



C. Heterogeneous Spiders 

Various spiders can be assembled experimentally 0], 
including those with distinguishable legs. Here we ana- 
lyze the coarse-grained properties of these 'lame' spiders. 

The bipedal lame spider is characterized by the max- 
imal separation s between the legs and by the hopping 
rates a and P of the legs, e.g., the a-leg hops to the right 
and left with the same rate a (whenever hopping is pos- 
sible) in the symmetric case. For the bipedal spider with 
s = 2, the diffusion coefficient is given by 



D2 



1 a(3 



2 a + P 



(19) 
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When a = 13 we recover the already known resuh telling 
us that the diffusion coefficient is 4 times smaller than the 
hopping rate. For a very lame spider (a ^ /3), Eq. (fT9|) 
gives D2 = a/2, so the diffusion coefficient is half the 
hopping rate of the very slow leg. 

To derive (fig]) , we first note that the probabilities (fTO]) 
satisfy 

dP 

-f =/3Q„ + aQ„_i-(a + /3)P„ (20a) 
at 



dQn 



= aP„+i + - (a + P)Qr. 



(20b) 



Here we have assumed that the left leg hops with rate a 
and the right leg hops with rate /3. (Recall that in one 
dimension, the order of the legs never changes.) 

Using Eqs. (|20ap - (j20bp we find that the mean-square 
displacement ([T^ evolves according to 



where 



dt 



E 



n— — 00 



a + /3 



Pn 



— {a — f3){u — v) 



(21) 



{n+l)Qn 



n— — 00 



From Eqs. ()20a|) - (|20bp we deduce that the quantity u — v 
obeys 

l{u-v) = ^~2{a + (3)iu~v) (22) 

Equation (I22|) shows that u — v quickly approaches to 
{a - /3)/[4(a + /?)]. Plugging this into ^ yields 



21 



a + P 



a + P 



which leads to the announced result, Eq. (fT9|) . 

For the biased bipedal lame spider, the drift velocity 
is given by a neat formula 



V2 = 



a/3 
a + P 



which resembles to p9|) . To establish 
the analog of Eqs. (|20ap - (l20ap . namely 



(23) 



one can use 



dPn 
dt 

dQr. 

dt 



= aQn-i - (3P„ 

= PPn - aQn 



(24a) 
(24b) 



Equations (|24a|) - (j24bp give the weights 



Further analysis of Eqs. (|24ap - (j24bp allows one to de- 
termine the diffusion coefficient 



1 ^ a'+P' 
2"^ (a + /3)3 



(25) 



We do not give a derivation of this formula since it can be 
extracted from earlier results by Fisher and Kolomeisky 
[l2i] who in turn used previous findings by Derrida jl3j. 

It is more difficult to compute the diffusion coefficient 
for the bipedal lame spider with maximal span s > 2. 
The results of Ref. [l^l do not cover the general case, al- 
though a proper extension of methods [HI, [l^ may solve 
the problem. For the symmetric bipedal lame spider with 
maximal span s > 2, we used an approach outlined in 
Appendix 1^ and obtained 



^P_ f _l 
a + P \ s 



(26) 



For s = 2, we recover equation ()T9|) . 

For the biased bipedal lame spider with maximal span 
s > 2 it is again simple to determine the drift velocity. 
Using an analog of p4ap - (|24bp one gets the weights and 
then the drift velocity is found from the relation 



2Vs = Pwi + {a + P) 2^ wi + aws 



s-l 

=2 



The outcome of this computation is 



Vs = aP 



a" 



(27) 



Specializing to a = /? = 1 (the I'Hospital rule allows to 
resolve an apparent singularity) one arrives at the ex- 
pression ^ for the velocity. 

Finally, the diffusion coefficient for the biased bipedal 
spider with arbitrary s is 



1 ^ - P 
— aP 

2 ^ a^~P' 



1 



s-l 



a 



s-l 



+ 



1 — s a 



a'' ~ P" a" - 

s+l _ ^sps 



(28) 



a* - P" a" - P" - 

This equation is derived in Appendix |Al 

Equation reduces to when s = 2; for s — ?> 
the diffusion coefficient can be re-written as 

1 (g + P){a' - a/3 + /3^)(a^ + 3q/3 + P^) 
'"2"^ (a2 + a/3 + /32)3 



Also when a = P — 1, equation (|28|) reduces to the 
expression Q for the diffusion coefficient. 



E 



W2 



n— — 00 



n— — oo 



a + P 

and relation V = {Pwi + aw2)/2 leads to 



E = 



a + P 



III. MULTI-PEDAL SPIDERS 

For the multi-pedal spider, L > 3, we must specify the 
constraint governing the separations between the legs. 
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FIG. 1: Illustration of a centipede, i.e. a spider with local 
constraint. All legs step to empty nearest neighbor sites at 
the same rate with adjacent legs staying within distance s 
from each other. 



A. Centipedes 

Here we consider centipedes or local spiders where the 
distance between the j^^ and (j + 1)*^* legs is at most s. 
See Fig. [T] for an illustration of such a centipede. In this 
case the total number of configurations is C = s^~^ since 
each of the (L — 1) spacings between adjacent legs can 
have s possible values. 



1. Main results 

Consider first spiders with s = 2. The configurations 
for the bipedal spider arc shown in the four possible 
configurations for the tripod are 



• o • o • 



(29) 



and generally there are 2^~^ possible configurations. 

Let D{L) be the diffusion coefficient of an L-leg spider. 
In the case of symmetric hopping (all rates are one) 



DiL) 



1 



4(L-1) 



(30) 



when s = 2. For L = 2, this of course agrees with our 
previous result: D{2) = D2 ~ 1/4. 

For the biased multi-pedal spider, the velocity is 



V{L) 



1 L+1 

2 2L- 1 



(31) 



The biased infinite-leg spider has a finite limiting speed! 
More precisely, V{oo) = 1/4, i.e., the infinite- leg spider 
drifts 4 times slower than the single-leg spider. The dif- 
fusion coefficient of the biased spider is 

3 (4L-3)![(L-1)!(£ + 1)!F 
~ 4 [(2i-l)!]3(2L + l)! ^^^^ 

Note that the diffusion coefficient of the infinite- leg spider 
vanishes. Asymptotically, 



D{L) 



3V27r 
128 



L 



(33) 



The above results dSHll-lISll) are valid when s = 2. We 
have not succeeded in computing V{L) and D{L) for arbi- 
trary L when the maximal separation exceeds two, s > 2. 

The velocity and the diffusion coefficient can be com- 
puted for centipedes with s > 2 when the number of legs 



is sufficiently small. The simplest quantity is the velocity 
of biased spiders. When s = 3, we computed the velocity 
V{L) of centipedes with up to seven legs: 

V{2) = 2/3 

V{3) = 26/45 w 0.5778 
1/(4) = 2306/4301 « 0.5362 
2257932864491452 ^ 
^ ' 4410656468591479 
V{6) w 0.4960476429 

V{7) w 0.4848259795 

(we have not displayed exact expressions for V{6) and 
V{7) which are the ratios of huge integers.) Note that 
for biased spiders with s = 3 one can guess the general 
expression pT|) from exact results for V{L) for small L; 
in contrast, no simple expression seems to exist for the 
velocity of biased spiders with s > 2. 

For symmetrically hopping spiders, we computed the 
diffusion coefficient when the number of legs is small. 
Here are the results for centipedes with s — 3 (the 
method used in calculations is described in Appendix E]) 



D{2) 
D{3) 

m 

D{5) 
D{6) 



1/3 

22/117 « 0.1880 
530/4059 « 0.1306 
145730406362990^ 
1457669284934841 

13157424169190800305558220463956878370565 
162454344889141072641777603974162004103911 
0.080991519 



In contrast to the neat formula (|30p characterizing the 
s = 2 case, the above numbers look intimidating. Fac- 
torizing the nominator and denominator of £'(6) reveals 
the presence of extraordinary huge factors and thereby 
excludes that it can be described by a formula like (|32p . 
let alone ([30)1 . Note that at least the asymptotic 
behavior predicted by Eq. (|30p remains valid for all s; 
for s = 3, in particular, we have D(L) ^ AL~^ with 
A sa 0.423 when L > 1. 



2. Mapping to the exclusion process for s = 2 

The derivations of above results are complicated since 
the number of configurations grows exponentially with L. 
Further, the transition rates are configuration dependent, 
e.g., for the four- leg spider configurations 

•••• • o • • • ••o«« 

evolve with rates 2, 3, 4 for symmetric hopping. (In 
contrast, for bipedal spiders the number of configurations 
grows linearly with s and the transition rates are simple.) 
All this makes the computation of the diffusion coefficient 
D{L) for arbitrary L very challenging. The pedestrian 
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calculation is feasible for small L, but even for L = 3, 
the framework based on rate equations like (|15ap ~ (|15cp 
is very cumbersome. 

Fortunately, spiders with local constraint and s = 2 
are related to simple exclusion processes (SEPs). This 
allows us to extract some predictions about spiders from 
previously known results about SEPs, and to employ the 
methods developed in the context of SEPs to situations 
natural in applications to spiders. 

We now demonstrate the remarkable connection be- 
tween centipedes with s — 2 and SEPs. As an example 
we show that the biased spider is isomorphic to the to- 
tally asymmetric simple exclusion process (TASEP) with 
open boundary conditions. To understand the isomor- 
phism, consider for concreteness the tripod. We can map 
configurations onto configurations 



00 



10 



01 



11 



(34) 



of the exclusion process on two sites with open boundary 
conditions. Here on site implies that there is no 
empty site between j**^ and (j -I- 1)'^* legs, while 1 implies 
that there is an empty site. A hop to the right of an 
internal leg in (^5)) corresponds to a hop to the left of a 
particle in Further, the hop of an extreme right leg 
corresponds to the addition of a particle to the extreme 
right position, and the hop of the extreme left leg corre- 
sponds to the removal of a particle from the extreme left 
position. The same mapping applies to any L. Thus in 
this TASEP each site i = 1, . . . , L — 1 can be occupied 
by a particle, and each particle hops to the left with rate 
one if this site is empty; further, a particle is added to 
site i = L — 1 with rate one if this site is empty, and a 
particle is removed from site 1 with rate one if this site is 
occupied. Thus we have shown that the s — 2 biased spi- 
der that moves to the right is equivalent to the TASEP 
with open boundaries in which particles hop from right 
to left. A similar mapping holds between the symmetric 
spider and the symmetric exclusion process. 

Derrida, Domany, and Mukamel fTi^ have shown that 
Eq. (pij) gives the flux in the TASEP; the isomorphism 
between the flux and velocity proves that the velocity of 
the biased spiders is given by ([51]) . This result was re- 
derived by other techniques, e.g., by a pure combinatorial 
approach [Tsj . The (much more c omp licated) derivation 
of the diffusion coefficient in Ref. [161 gives 



First we recall that the symmetric spider with L legs 
hopping in both directions with rates equal to 1 is equiv- 
alent to a symmetric exclusion process on L — 1 sites with 
open boundary conditions. For this SEP, all rates (i.e. 
hopping rates in the bulk, entrance and exit rates at the 
boundaries) are equal to 1 . This Markov process satisfies 
detailed balance and is at equilibrium; in particular, the 
mean current, i.e., the velocity of the spider, vanishes 
identically. The variance of the current corresponds to 
the diffusion constant of the spider. This variance can be 
calculated as follows. 

Consider now a symmetric exclusion process of length 
L — \ with open boundaries and arbitrary addition and 
removal rates at the boundaries. The system is driven 
out of equilibrium by particles entering and leaving at 
the boundaries. In the bulk, each particle hops with rates 
1 to the right and to the left (if the corresponding sites 
are empty); a particle enters at site 1 with rate a and 
leaves this site with rate 7; similarly a particle enters 
another boundary site L — 1 with rate 5 and leaves this 
site with rate /3. Generically, these unequal rates lead to 
a current. The mean value of this current is given by (see 
e.g. Ref. [13) 




0+S 



Q+7 



L 



1 

a+7 



1 



(35) 



The equilibrium conditions correspond to a = ^ — [3 — 
5 = 1 and J = 0. We now choose the rates on site 1 
as follows a = exp(|) and 7 = exp(— |) and we keep 
(3 = 5 = 1 &t site L — 1. Then, the current is given by 



J 



tanh ; 



2L-3 



(36) 



The Markov matrix of this process satisfies the gener- 
alized detailed balance condition given by equation (|B2p 
of the Appendix [B] with = ±1 if a particle enters at 
site 1, or exits from site 1 (y = otherwise). We can 
then use the fluctuation-dissipation formula (see (IB4|) in 
the Appendix |B]) wich tells us that the fluctuation of the 
current at the first site is given by 



D 



dJ_ 

a7 



1 



e=0 4(L-1) 



(37) 



in accordance with equation ([30 



3. Derivation of ([50)) 

For the symmetric spider, it should be possible to com- 
pute the diffusion coefficient ([50[) by using the technique 
of Ref. [1^. This technique (based on an extension of 
a matrix technique) is very advanced. The final result 
([30[) looks much simpler than its biased counterpart ([32[) . 
Hence we have sought another derivation, and we have 
found an intriguingly simple proof of (j30D based on the 
fluctuation-dissipation formula (see Appendix [B[) . 



4. Mean-field approximation for s > 3 

Simple exclusion processes have been thoroughly inves- 
tigated (see books and reviews Jl8, 19, 20, 21]). Hence 
one can extract the results about spiders from already 
known results about SEP. Unfortunately, for spiders with 
local constraint the mapping onto SEP applies only when 
s = 2. The spider with L legs and arbitrary s can be 
mapped onto an exclusion-like process with L — 1 sites 
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and with open boundaries. In this process the maximal 
occupancy is limited, namely the number of particles in 
each site cannot exceed s — 1. The dynamics is simple: 
one chooses sites with rate one and moves a particle to 
the site on its left; nothing happens if the chosen site 
was empty or the site on the left was fully occupied. One 
also adds particles to site i — L — 1 and removes from 
site i = 1, both these processes occur with rate one; the 
addition is possible as long as site i = L — 1 is not fully 
occupied (contains no more than s — 1 particles) . Unfor- 
tunately, this neat process has not been solved exactly 
but it can be studied by a mean-field analysis. 

To simplify the analysis, we consider centipedes with 
infinitely many legs. We assume that the distance be- 
tween adjacent legs cannot exceed s. We further assume 
that the spider's motion is biased, and limit ourselves to 
a (mean-field) computation of its velocity V^-^^ . 

First we map the spider onto the generalized asym- 
metric exclusion process with at most s — 1 particles per 
site. We then write Xj for the density of sites with j 
particles; this is just the density of gaps of length j -I- 1 
between adjacent legs of the spider. The possible values 
are j = 0,...,s — 1. Writing the evolution equation for 
Xj and setting Xj = we obtain 

{xj-i - Xj){l - xq) - {xj - x.j+i){l - Xs-i) = (38) 

when 1 < J < s — 2. Similarly from xq — and i^-i = 
we get 

xi(l - a-s_i) - xo(l - Xq) = (39a) 
a;s_2(l - a;o) - a;s-i(l - a;s-i) = (39b) 

The obvious normalization requirement is 



3=0 



(40) 



As a warm up, consider the first non-trivial case s = 3. 
Due to normalization, it is sufficient to use (|39a|) - (j39bp . 
Writing xq = x and X2 = z, we have xi = 1 — x — z from 
(HO)) , and (pga|l - (l39bl) become 

{1-x- z){l- z) = x(l-x) (41a) 
(1 - .T - z)(l - x) = z(l - z) (41b) 

These equations are actually identical; solving any of 
them we arrive at 



^ z+v/z(4-3z) 
2 



(42) 



To compute velocity we return to the original formu- 
lation. A leg of the spider moves with rate 1 if the site 
ahead is empty and if the leg behind is one or two steps 
behind. The former event happens with probability 1 — x, 
while the latter occurs with probability I — z. Thus the 
velocity is 



y(3) = {l-x){l-z) 



Using (|42l) we get 



V 



-(1-z) [z + ^z{A-3z) 



(44) 



We should select maximal velocity. The maximum of 
V{z) given by (|44|) is reached at z = 1/3, and it reads 



1/(3) 



(45) 



At the state corresponding to the actual (maximal) ve- 
locity all densities are equal: xq = xi = X2 = 1/3. 

The situation for s > 3 is also simple. Analyzing re- 
currence ([55]) one finds that for all < j < s — 1 the 
solution is a shifted geometric progression 



Xj = A 



A 



Xq 



1 - X. 



(46) 



Plugging (|46p into Eqs. (|39ap - (j39bp one achieves the con- 
sistency if either A = or A = 1. In the latter case the 
densities are the same, and hence they are all equal to 
due to normalization requirement (|40p . The straight- 
forward generalization of (|43|) is 



1/(^)-(1-xo)(1-:e,-i) 



and therefore 



(47) 



(48) 



In the complimentary case of A = the analysis is a 
bit more lengthy. However, the final result is the same. 
Here is the proof. Since Xj = BX^ , equation gives 
A = {1 — B)/{1 — _BA*^^), which can be re-written as 



B = 



1 - A 
1-A^ 



Further, (|T7|) becomes 

V^(") = {I- B){1- BX'-^) 
Using (|49|) we recast ([50]) into 

V s-l \ 2 



y(^)(A) - A 



1-A" 
l-A'* 



(49) 



(50) 



(51) 



(43) 



The maximum of V'-'\X) is achieved at A = 1. Thus the 
velocity is indeed given by Eq. (|48)) . 

The above elementary analysis is mean-field as we have 
assumed the validity of the factorization. The answer is 
trivially exact for s = 1, and it is known to be exact for 
s = 2. For s = 3, we calculated velocities exactly for 
small centipedes, see Sect. IIII A II The limiting L ^ oo 
value obtained from simulations T^^^) ~ 0.4189 is close 
to the predicted mean- field value y(3) ^ 4/9 ^ 0.4444. 
Overall, the assumed factorization is not exact when s > 
3. Note that a model which differs from our model only 
in the hopping rules has been solved exactly [l^l, but 
there the stationary state is a product measure. 



5. Lame spiders 

Finally we investigate lame centipede spiders whose 
extreme left leg hops to the right with rate a and the 
extreme right leg hops to the right with rate (3. The above 
mean-field analysis shows that the velocity of the extreme 
left leg \s a{l — xq) and the velocity of the extreme right 
leg is /3(1 — Xs-i)- As long as these velocities exceed the 
bulk velocity (|48|) . the actual gap density xq at the left 
end and Xg^i at the right end will be higher than their 
bulk values, so the spider will move with velocity 
This occurs as long as a(l — s^^) and /3(1 — s^^) exceed 
(1 — s~^)^, i.e. a,P > 1 — s^^. When at least one of 
the rates is smaller than the threshold value, different 
behaviors emerge. Overall, the speed of the infinite-leg 
spider exhibits an amusing dependence on the rates a 
and /3: 

( (1 - s^^ f for a,/3 > 1 - s-^ 
V^"^ ^ \ Ws{a) for a < /3, a < 1 - (52) 
[ W,{f3) for f3<a, (3 < I - s"^ 

Thus if at least one of the two extreme legs has the in- 
trinsic speed less than 1 — s~^, the speed of the entire 
spider is solely determined by the slowest leg. 

To determine Ws{(3) we note that velocity on the right 
boundary is 

V^Pil- xs-i) = fi{l - B\'-') (3 (53) 

where in the last step we have used (|49| . Equating the 
velocity given by Eq. ([53]) with the velocity in the bulk 
given by Eq. (|5ip we find 

/? = A^-^ (54) 

Thus the velocity is given by ((53)) or ([5T|) . where param- 
eters are connected via ([54]). 

Explicit results can be obtained for s up to s = 5. For 
s = 2 we recover the celebrated result 

W2W)^[3{l-P) (55) 

For s — i the final expression is still compact 



T^3(/3)=/3(l-/3)ii^^^^ (56) 

with b — j3/{l — (3). For s = 4 the result is quite cum- 
bersome 

Wm = ^\. A=-A--A-i-- (57) 
where we have used the shorthand notation 

A = (28 + 1086 + 12v/9 + 426 + 8152}^^^ 
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FIG. 2: Illustration of a spider with global constraint. The 
legs can step independently to nearest neighbor empty sites 
within a distance S from each other. 



B. Global Constraint 

For spiders with L legs and maximal span S between 
any two legs (see Fig. [2]), the global constraint rule hmits 
the maximal distance between the extreme legs and the 
exclusion condition implies that S > L — 1. A spider 
with maximal distance S" = L — 1 is immobile, so we shall 
tacitly assume that S > L. It is also useful to keep in 
mind that for a spider satisfying the local constraint rule 
the maximal span is (L — l)s if the maximal distance 
between the adjacent legs is s; for the bipedal spider 
S = s. A spider with global constraint is equivalent to the 
exclusion process on a ring, where each leg is interpreted 
as a particle and the total number of sites is equal to 5-1-1. 
For such a process with periodic boundary conditions, a 
key property of the stationary state, which holds both 
in symmetric and biased cases, is that all configurations 
have equal weight (Toj . 



1. Configurations 

To count the total number of configurations, we set, 
as usual, the origin at the position of the extreme left 
leg, see e.g., (p9t : this allows us to avoid multiple count- 
ing of configurations which differ merely by translation. 
We then note that the other L — 1 legs can occupy sites 
1, . . . , S. Thus the total number of configurations is 

e(i,^)-(^!^) (58) 

In the stationary state, the weight of a configuration is 
thus given by w = 1/6. 

Let us now calculate the total number INf(L, S) of "o 
pairs in all configurations. Each configuration begins 
with a string 

.o_^. (59) 

a 

where a = 0,1, . . . , S — L + 1. Disregarding the part up to 
the second leg maps configurations of the type ([55]) with 
fixed a to configurations of the spider with L — 1 legs and 
maximal span S~a. The total number of "o pairs in these 
latter configurations is 'N{L — 1, S — a). Configurations 
of the type ([55]) have of course an additional •o pair at 
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the beginning (when a > 0). Therefore 

S-L+l S-L+1 

'N{L,S)= ^{L-l,S-a)+ e{L-l,S~a) 

The latter sum is simplified by using (|58p and the identity 



qj + 1 



Thus we arrive at the recurrence 

S-L + l 

J<{L,S)= Y ^{L-l,S~a) 

a=0 



S-1 
L-2 



(60) 



The solution (found by a generating function technique 
or verified by mathematical induction) reads 



N(L,5) = L 



S-1 
L - 1 



(61) 



Since all configurations have equal weight in the sta- 
tionary state, the velocity of the biased spider can be 
expressed by the total number [NT of "o pairs as 



V = L 



-1^ 



which then leads to 



V = 1- 



L- 1 



(63) 



(Note that the velocity is zero in the unbiased case.) It is 
more involved to calculate the diffusion coefficient, which 
we obtain below separately for the unbiased and for the 
biased case. 



2. Symmetric hopping 

The diffusion coefficient of the symmetric spider is 
given by 



(64) 



This equation is obtained by applying the fluctuation- 
dissipation relation, which is valid because the dynam- 
ics of the symmetric spider statisfies detailed balance (in 
other words, the symmetric spider is a system in thermo- 
dynamic equilibrium) . The Einstein relation then implies 
that D (X V, where the velocity V of the biased spider 
is given by the expression (|62p . The extra factor 
between ([M)) and comes from the fact that D is the 
diffusion coefficient of the left leg (or of the spider's center 
of mass) whereas V is the mobility of the biased spider 
where all the legs are asymmetric. Using Eqs. ([58| and 
(pT|) we recast ^ to 



L 



1 - 



S 



(65) 



For the bipedal spider there is no difference between local 
and global constraints. Using L — 2 and S* = s we find 
that Eq. ([GS]) indeed turns into Eq. Q. The calculation 
of D presented here is self-contained; we notice that the 
expression ([65]) can also be found as a special case of a 
general formula derived in [23| for the diffusion constant 
of a partially asymmetric exclusion process. 



3. Biased hopping 

While the velocity ([55)1 is easily computable, the dif- 
fusion coefficient was calculated in [2^ using a matrix 
Ansatz (see [l^ and [1^ for a more general formula). 
The result is 



D{L,S) 



2(25 -2L+1) \2L-lJ \L-1 



2S -I 



(66) 



For a given number of legs, the diffusion coefficient of 
the most clumsy spider is 



(67) 



(62) while the diffusion coefficient of the most agile spider is 



D{L, oo) 



-)2L~2 



2L 



L \ L 

When L ^ 1, the diffusion coefficient (pS)) scales as 



(68) 



(69) 



More generally, the diffusion coefficient D{L, S) also de- 
creases as (7r/16L)^/^ when 1 < L < S^^^. 



C. Heterogeneous Spiders 

Each leg of a heterogeneous (lame) spider may have its 
own hopping rate. The bipedal lame spider was studied 
in Sec. Ill CI One can find explicit expressions for the ve- 
locity and the diffusion coefficient of the lame tripod and 
perhaps for the lame spider with four legs; the general 
solution for an arbitrary L is unknown. 

Lame spiders are tractable if only one or two legs have 
different hopping rates. Below we consider lame spiders 
whose extreme legs are affected. For concreteness, we 
focus on spiders with local constraint and s = 2. The 
analogy with the TASEP with open boundary conditions 
still applies, the only modification is that the particle is 
removed from site 1 with rate a and the particle is added 
to site L — 1 with rate (3. The flux in such a system was 
found in Ref. [l^l; this gives us 



CL_i(a,/3) 



(70) 
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where we used the shorthand notation 

N 



p{2N-l-py. a-P~^ - p-P~^ 



p 



^ Nl{N-p)l a-i-/3-i 



Plugging Cq — 1 and Ci = a^^+P^^ into ([70)1 we recover 
the expression (|23p for the velocity of the bipedal lame 
spider; the velocity of the lame tripod is 



^3 («,/?) 



The speed of the infinite-leg spider exhibits an amusing 
dependence on the rates a and f3: 




for a > 1/2, P > 1/2 
a) for a < P, a < 1/2 
/3) for P <a, /3 < 1/2 



(71) 



Thus if both rates exceed 1/2, the speed attains a univer- 
sal (independent of the rates) maximal value Voo = 1/4. 
On the other hand, if at least one of the two extreme 
legs has the intrinsic speed less than 1/2, the speed of 
the entire spider is solely determined by the slowest leg. 

A general explicit expression for the diffusion coeffi- 
cient is unknown. There are two special cases, however, 
in which the diffusion coefficient was explicitly calculated 
[l6| . One is the homogeneous spider (a = /3 = 1) when 
D{L) is given by Eq. ([5^ : another particular case corre- 
sponds to a -|- /3 = 1 when the diffusion coefficient is 



Dr 



L-2 

E 

k=0 



2(2fc)! 
kl{k + l) 



fe+i 



(72) 



with Voo given by Eq. ([TTjl ; since ([72|) is valid on the line 
a + P = 1, we have Voo = a(l — a) = P{1 — P). As a 
consistency check one can verify that equations (j72p and 
(|25| do agree : setting L = 2 in the former and a + P = 1 
in the latter we indeed obtain the same result. 

The behavior of Dl for the spider with many legs is 
again amusing. For the infinite-leg spider, Eq. ([7^ yields 



Doo = ^aP\a 



P\ when a + P = 1 



(73) 



Thus on the line a + P = I, the diffusion coefhcient van- 
ishes only when a — P — 1/2. 

The behavior of the diffusion cocfhcicnt for the infinite- 
leg spider is particularly neat, and it had actually been 
understood (in the context of the TASEP) for arbitrary 
a and p. Derrida, Evans and Mallick [l^ found that 



Do 
14c 





(1 - 2a)/2 
(1 - 2P)/2 
(1 - 2/3)/3 



for a > 1/2, P > 1/2 

for a < P, a < 1/2 

for /3 < a, P < 1/2 
for P < 1/2 



(74) 



The discontinuity on the symmetry line a = P < 1/2 is 
especially striking. 



IV. QUICK SPIDERS 

In the previous sections we have considered the sim- 
plest possible gaits when the spider's legs can step only 
to the neighboring sites. In this section we briefly ex- 
plore the behavior of quick spiders. These spiders (in- 
troduced in Ref. [l3|) differ from previously discussed 
spiders, namely the legs of a quick spider can jump over 
several lattice sites at once. The only requirement is to 
stay within distance S from the other legs. Hence quick 
spiders can be in the same states as the corresponding 
global spiders, but more transitions are possible between 
the states of the quick ones. 

The simplest quick spider has two legs always next to 
each other (L = 2, S" = 1). Although such a global spider 
cannot move, a quick spider can put one leg ahead of the 
other and can walk this way. Its motion is completely 
equivalent to a simple random walk, hence its diffusion 
coefficient is D — 1. This is generally true for quick 
spiders with L legs and maximal distance S = L — 1 



D{L,S = L-1) = 1 



(75) 



We also computed the diffusion coefficient of bipedal 
quick spiders with arbitrary S. We found 



D{2,S) = 



S{S+1){2S+1) 
6 



(76) 



This expression can be derived using the general formula 
([64| . For the bipedal spider we can label various config- 
uration by the distance 1 < £ < S between the legs 

(77) 



Take the left leg. It can jump to the left up to distance 
S — £; the corresponding displacements of the center of 
mass are Aa; — —i/2 with 1 < i < S — £. The left leg can 
also jump to the right. The displacements are Ax — i/2 
with 1 < i < i — 1, and once it overtakes the right leg, 
Aa; — (£ + i)/2 with 1 < i < S. Taking also into account 
that all weights are equal, wg — 1/5, and recalling that 
jumping of the right leg will give the same contribution, 
we recast (l64l) into 



1 ^ 

e=i 



4=1 



i=l 



4=1 



(78) 



Computing the sum yields the announced result (I76p . 



V. INTERACTING SPIDERS 

In experiments Q, thousands of spiders are released, 
yet their density is usually small. Naively, one can an- 
ticipate that spiders are essentially non-interacting. This 
is correct in the earlier stage, t < t^, but eventually spi- 
ders "realize" the presence of other spiders, and their 
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behavior undergoes a drastic change from diffusive to a 
sub-diffusive one. This intermediate stage proceeds up to 
time t* when spiders explore the entire system and then 
the diffusive behavior is restored, albeit with a smaller 
diffusion coefhcient 2). Here we compute D and estimate 
the crossover times and t* . 

Let N spiders be placed on the ring of size S. We 
assume that the spider density n — N/S is low, n <C 1; 
equivalently the typical distance (n~^ lattice spacings) 
between neighboring spiders is large. 

Imagine that we know the diffusion coefficient D of an 
individual spider (e.g., for the bipedal spider with s — 3 
we found D — 1/3 when each leg hops symmetrically 
with rates equal to one) . Each spider covers around ^/Dt 
lattice sites, and equating \Jlytl = n"^ we arrive at the 
estimate of the lower crossover time 



The behavior is sub-diffusive in the intermediate time 
range, <C t <C t*. It is characterized by the {Xt^l^ 
growth of the covered line [1^; this so-called single- file 
diffusion has numerous applications [l^, • The ampli- 
tude A is found by matching [Dt^Y = (At,)^/'* which 
in conjunction with ([79|) yield A = D/n^. 

The final behavior is again diffusive. In the long time 
regime, t ^ t* , we may interpret each spider as an ef- 
fective particle hopping to the right or left with rates D. 
The interaction between spiders is essentially equivalent 
to exclusion interaction between particles, and hence the 
system reduces to the SEP. We can therefore use (p5)) 
where we should replace L by TV, and we must also mul- 
tiply the result by D since spiders effectively hop with 
rates D rather than one. The term in the brackets in 
Eq. (j65p reduces to 1 — n; we can replace it by one since 
n <C 1. Therefore Eq. ([65)1 becomes 

D = N-^D (80) 

Thus exclusion interaction greatly reduces the diffusion 
coefficient. This strong cooperative effect emerges even 
when the density is arbitrarily small, the only require- 
ment is that there are many spiders, ^ 1. 

The upper crossover time t* is found by equating 
P^*)i/2 ^ (Ar)i/4. We arrive at 

r = ^ = Nh, (81) 

Thus the analogy with SEP essentially solves the prob- 
lem in the practically important limit when the spider 
concentration is low. Neither memory nor the gait play 
any role, one must merely use the diffusion coefficient D 
corresponding to the actual gait and computed under the 
assumption that the lattice sites are in the product state. 
One should remember, of course, that the SEP regime is 
achieved when t > t*; at much earlier times t > t<t, the 
spiders mostly hop on the product, and therefore the as- 
sumption of full attachment can become problematic. 



VI. DISCUSSION 

A single spider is a self-interacting object. There are 
two sources of interaction between the legs: (i) exclusion 
(no more than one leg per site), and (ii) legs cannot be 
too far apart. Is it possible to represent a spider as an 
effective single particle? The answer is yes — at least in 
simple situations, one can treat a spider as a diffusing 
particle. It is far from trivial, of course, to compute the 
diffusion coefficient of this particle. Fortunately, natural 
models of spiders are related to simple exclusion pro- 
cesses. In the course of this work we had an advantage of 
utilizing some beautiful results and powerful techniques 
developed in the studies of simple exclusion processes. 

Our models certainly do not take into account all the 
details of an experimental situation Q. For instance, we 
assumed that the re-attachment of a leg is very quick, so 
the process is controlled by detachment. Hence spiders 
remain fully attached and never leave the surface. This 
assumption is important as our analysis has relied on the 
permanent presence of spiders on the surface. Relaxing 
this assumption does not make the problem intractable 
— indeed in recent analyzes of molecular motors the com- 
plete detachment (unbinding) from cytoskeletal filaments 
is allowed, see e.g. j29, 30, 31, 32, 33]. Further, our anal- 
ysis of the many-spider situation in Sec. |V]treats the low 
density case; the analogy with SEP allowed us to handle 
the problem but the assumed permanent presence of the 
spiders is particularly questionable in this case. 

Perhaps the most serious limitation of our analysis is 
the disregard of memory — in experimental realizations 
spiders often affect the environment which in turn 
affect their motion. The non-Markovian nature of this 
problem calls for a set of new techniques even in the case 
of a single spider. In one dimension, the influence of 
memorycan be probed analytically for a single bipedal 
spider @], and the replacement of a self- interacting spi- 
der by an effective particle remains valid, though this 
effective particle becomes an excited random walk which 
distinguishes visited and unvisited sites. 

Finally we note that the SEP and its generalizations 
occur in various biolo gica l problems ranging from motion 
of molecular motors pol . [sol . [3ll . [33 , |33| | to protein syn- 
thesis [33, [SI, [3I, [s^l ■ Some models of protein synthesis 
resemble complicated models of spiders. Another intrigu- 
ing connection is between spiders and cooperative cargo 
transport by several molecular motors [ssj . 
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APPENDIX A: MASTER EQUATION AND 
FLUCTUATIONS 

In this Appendix, we explain the general formalism, 
inspired by Ref. [l^, that allows one to calculate veloc- 
ities and diffusion contants, and we use this method to 
derive equation ([28|) . 

A spider can be viewed as a homogeneous Markov pro- 
cess with a finite number of internal states. The dynam- 
ics of the spider is encoded in a Markov Matrix M, where 
the non-diagonal matrix element M{C, C) represents the 
rate of evolution from a configuration C" to a different 
configuration C. The quantity — M(C, C) is the exit-rate 
from configuration C. The master equation for Pt{C), 
the probability of being in configuration C at time t, is 
then given by 



-Pt{C) = Y,M{C,C')Pt{C'). 



(Al) 



We now define Ij as the absolute position of the spider's 
left leg, knowing that at time t = Yt — Q. Between 
t and t + dt, Yt varies by the discrete amount -1-1,0 or 
— 1 that depends on the configuration C at t and on the 
configuration C a.t t + dt. The Markov Matrix M can 
then be decomposed in three parts corresponding to the 
three possible evolutions of Yt : 

M{C,C') = Mo{C,C') + Mi{C,C')+M^i{C\C') . (A2) 

For example, Mi{C,C') represents the transition rate 
from a configuration C" to C with the left leg moving 
one step forward (this matrix element vanishes other- 
wise); M-i corresponds to transitions for which the left 
leg moves one step backwards; Mq encodes transitions 
in which the left leg stays still. We call Pt (C, Y) the 
joint probability of being at time t in the configuration 
C and having Yt = Y. A master equation, analogous to 
equation (|Aip . can be written for Pt{C, Y) as follows : 



dt 



PtiC,Y) = J2{Mo{C,C')Pt{C\Y) 
c 

+ M^{C,C')PtiC',Y-l) 

+ Af„i(c,c')P*(c',y + i)) 



(A3) 



In terms of the generating function Ft{C) defined as 

oc 

Ft{C)^ J2 e^^Pi(C,r), (A4) 

Y=-oo 

the master equation IjASp takes the simpler form : 

^F,(C)=^M(A;C,C")F,(C"), (A5) 

c 

where M (A; C, C"), which governs the evolution of Ft{C), 
is given by 



We emphasize that M{\), is not a Markov matrix for 
A ^ (the sum of the elements in a given column does 
not vanish). 

In the long time limit, t — > oo, the behaviour of Ft{C) 
is dominated by the largest eigenvalue /i(A) of the matrix 
A/(A). We thus have, when t — > oo, 



[e^y^)=Y^Ft{C)^e 



(A7) 



This result can be restated more precisely as follows 



lim - log( e^^* ) = m(A) 

t^oo t 



(A8) 



The function /i(A) contains the complete information 
about the cumulants of Yt in the long time limit. For 
example, the velocity V and the diffusion coefficient D 
of the spider are given by 



V 



D = hm 

t—*oo 



lim<M = M^ 

t^oo t dA 



m"(o) 



2t 



(A9) 
(AlO) 



One therefore needs to calculate the function /i(A). For 
simple problems such as the bipedal spider with s = 2, 
/i(A) can be determined exphcitely (because M(A) is a 
2 by 2 matrix). In general, the most efficient technique 
is to perform a perturbative calculation of /i(A) in the 
vicinity of A = (recall that /i(A) vanishes at A = 0). 
This perturbative approach is very similar to the one 
used in Quantum Mechanics, the major difference being 
that M{X) which plays the role of the Hamiltonian is not, 
in general, a symmetric matrix and its right eigenvectors 
are different from its left eigenvectors. By definition, we 
have 



m(a)|m(a)) = 

(MA)|Af(A) = MA)(MA)|. 



(All) 



M(A) = Afo + e^Afi -1- e~^A//_i 



(A6) 



Using equations (IA6p , (jA9p , and (jAlOp , we can write the 
following perturbative expansions in the vicinity of A = 0, 

M(A) = Af + A(Af,-Af_,) + ^!i^^il±^... 

/i(A) = V\ + D\^ + ... 
|/i(A)) = |0)+A|1)+A2|2) + ..., 
(MA)| = (0| + A(1| + A2(2| + ... 

where M is the original Markov matrix of the system, | ) 
is the stationary state and ( | = (1, 1, . . . , 1) is the left 
ground state of M . We now substitute these perturbative 
expansions in (|A11[) and identify the terms with the same 
power of A. Using the left eigenvector (/i(A) |, we obtain 

(0|A/ = 0, (A12) 
(l|Af = y(0|- (0|(Afi-Af_i), (A13) 

(2|Af = i?(0|-i(0|(A/i-|-A//_i) 

-^F(l|-(l|(Afi-Af_i). (A14) 
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Multiplying these equations by the right ground state | ) 
of M, and using the fact that M| ) = and ( | ) = 
1, the following formulae for V and D are derived as 
solvability conditions for Eqs. (|A12p - (|A14|) : 



V 
D 



;0|A/i-Af_i|0) 
;i|Mi-M_i|0> 



(A15) 



1 



+ -(0|Afi + Af_i|0) - V(1|0). (A16) 



We observe that in order to calculate V we only need to 
know the ground state of M. However, the calculation 
of D requires the knowledge of ( 1 1 , obtained by solving 
the linear equation (|A13p . We remark that similar ex- 
pressions can be obtained starting from the expansion of 
right eigenvector | ^{X) ). 

We now specialize this framework to the case of the 
heterogeneous bipedal lame spider with s internal states. 
The Markov Matrix is then an s x s matrix M = Mq + AIi 
since M_i vanishes identically. The matrix Mq is given 
by 



f-a 

a -{a + P) 

(3 -{a + 13) 







/3 -(a + /3) 

/3 -aj 



and the matrix Mi is 



Ml 



V 

The stationary state of M is | 

a — 13 



1 
0/ 



Pk 



(3^ 



{POiPi, ■ ■ ■ ,Ps-i) with 
s-k-ipk fQj, fc = o,...,s-l. (A17) 



This expression, together with (|A15[) . leads to the for- 
mula l|27p for the spider velocity. 

In order to derive the expression of the diffusion co- 
efficient, we need to solve equation (jA13l) . One can 
verify that the solution of this equation is given by 
( 1 I = {qo,qi, ■ ■ • where 



gfc = (fc + 1) 



V — a al3 — aV 
(3-a ^ {/3 - a)2 



(A18) 



for fc = 0, . . . , s — 1 . Inserting equations (|A17p and (jA18|) 
into the general expression (jA16p leads to the for- 
mula ([28]) . 

We also used the above method to determine the veloc- 
ity and the diffusion coefficient for centipedes with s = 3. 
The results (Sect. IIIIA1|) were obtained by explicitly 
constructing the matrices Mq,Mi, and M-i, and per- 
forming exact computations using Maple. These compu- 
tations are feasible when the number of legs is sufhciently 
small. (The total number of configurations is 3^~^ for 
centipedes with s = 3, and hence the order of matrices 
Mq, Ml, M_i quickly grows with L.) 



APPENDIX B: GENERALIZED DETAILED 
BALANCE RELATION 

For the symmetric spider, the three matrices A/q, Mi 
and Af_i, introduced in (jA2[) to take into account the 
total displacement of the spider, satisfy the following de- 
tailed balance relation 



MyiC C')P'"^iC') = Af_y(C", C)P'"i{C) 



(Bl) 



where the equilibrium measure is denoted by P'^'^ and 
y = 0,±1. Equation (|Bip implies that the velocity of 
the spider vanishes. Consider now a spider driven out of 
equilibrium with a non- vanishing mean velocity. Suppose 
however, that for the model under consideration there 
exists a real number e such that the following generalized 
detailed balance relation is satisfied 

Afy(e, e')P'^'(e') = M-v{&, e)F^«(e) exp(ey) , (B2) 

Here again, P'^'^ is the equilibrium measure corresponding 
to the symmetric spider. From relation (jB2p it is a matter 
of elementary algebra to prove that the spectra of Af (A) 
and of M{—e — A) are identical. Therefore 



A). 



(B3) 



This relation, which is a special case of the general 
Fluctuation Theorem valid for systems far from equilib- 
rium [3^ liO, Hlj , was derived for stochastic systems by 
Lebowitz and Spohn (4^ . Close to equilibrium, when 
e ^ 1, we can expand equation (|B3P for small A and e. 
We find at lowest order 



,"(0) = ^ ^.e. 



e 



(B4) 



which is nothing but the classical fluctuation-dissipation 
relation between diffusion and mobility. 
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